Fluctuations in the level density of a Fermi gas 
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We present a theory that accurately describes the counting of excited states of a noninteracting 
fermionic gas. At high excitation energies the results reproduce Bethe's theory. At low energies 
oscillatory corrections to the many-body density of states, related to shell effects, are obtained. 
The fluctuations depend non-trivially on energy and particle number. Universality and connections 
with Poisson statistics and random matrix theory are established for regular and chaotic single- 
particle motion. 
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The level density is a characteristic property of ev- 
ery many-body quantum mechanical system. Its precise 
determination is often a key ingredient in the calcula- 
tion of different processes, like compound nuclear decay 
rates, yields of evaporation residues to populate exotic 
nuclei, or thermonuclear rates in astrophysical processes. 
The first and main step towards the understanding of 
the density was given by H. Bethe more than sixty years 
ago. He considered A noninteracting fcrmions moving in 
a mean-field potential. By this model he showed that 
the number of excited states of the many-body (MB) 
system contained in a small energy window dQ at energy 
Q, p MB (A,Q)dQ, is 



p MB (A,Q) 



48 Q 



exp \/-7r 2 p Q . 



(1) 



Here ~p is the average single-particle (SP) density of states 
at Fermi energy e F , and Q is measured with respect to 
the ground state energy. For a given potential, p is in 
general a function of A. This expression assumes that 
the excitation energy is, on the one hand, large compared 
to the SP mean level spacing 8 — ~fT x and, on the other, 
small compared to e F (degenerate gas approximation). 
Most of the present knowledge concerning the density of 
states is based on refinements of Bethe's result. 

In Eq. JTJ all the information concerning the SP spec- 
trum is encoded in a unique parameter, p, that describes 
its average behavior. The exact dependence of p MB on A 
and Q is sensitive, however, to the detailed arrangement 
of the SP energy levels around the Fermi energy. Af- 
ter Bethe's work, more accurate calculations of p MB were 
made Schematic shell corrections related to a peri- 
odic fluctuation of the SP density were computed in Q 
(see also Q), introducing the so-called backshifted Bethe 
formula. Several modifications of Eq. have been done 
to match the experimental results. These models take 
into account, for example, shell effects, pairing correc- 
tions and residual interactions |5|, introducing a multi- 
tude of coexisting phenomenological parameterizations. 
The present status of the understanding do not allow to 



draw a clear theoretical picture of the functional depen- 
dence of p M B with A and Q. 

Within an independent particle model we demonstrate, 
through a general theory, the role played by shell effects 
in the MB density of states. Superimposed to the smooth 
growth of the density described by Eq. Q , there are os- 
cillatory corrections. To lowest order, these corrections 
arc directly related to energy fluctuations of the system, 
that may be expressed in terms of sums over the clas- 
sical periodic orbits. We make a statistical analysis of 
the density fluctuations, in particular of their probability 
distribution, typical size, and temperature dependence 
0. We show that the amplitude of the shell corrections 
grows linearly at low temperatures, and decays as T~ x 
when T 3> E c /2ir 2 , where E c is the energy conjugate 
to the time of flight across the system (temperature is 
measured in units of Boltzmann constant). The statis- 
tical properties of the density fluctuations strongly de- 
pend on the regular or chaotic nature of the underlying 
classical SP motion. Universal properties for each class 
(regular or chaotic) are found for T <C E c /2ir 2 . Besides 
the fluctuations, we find and discuss corrections to the 
smooth part of the density in the exponent of Eq. . The 
present analysis considers general properties of fermionic 
systems. Applications to the specific problem of the nu- 
clear density and to the interpretation of the experimen- 
tal results will be done in a separate contribution. 

The MB density of states is defined as 



Pmb (a, e) = S( - A - *o 6 ( E - E 



(2) 



N j 



where E N<j = J2i 
figuration oi N = 



is the energy of the j-th SP con- 
Si n j,i particles, ej the SP energies, 
and rij,i = 0, f are the occupation numbers. 

When the SP spectrum consists of equidistant levels 
separated by 8, e.g. the spectrum of a one dimensional 
harmonic oscillator, an exact answer to the counting 
problem exists. The MB excitation energies are integer 
multiples of 8, each level having a nontrivial degener- 
acy. It is easy to see that, for a given level defined by 



an integer, the degeneracy is the same as the number of 
ways into which the integer can be decomposed as a sum 
of integers (the partition number). An asymptotic ap- 
proximation to this well-known mathematical problem 
was obtained by Hardy and Ramanujan, and later on 
Rademacher found a convergent series Q. Expressing it 
in terms of an expansion in terms of pQ, the result reads 

log (A%m = \J~**PQ - log pQ) - (3) 

^<w-' -(51- a) w- 1 +»«»>-•> 

The first two terms of this expansion reproduce Bethe's 
formula. 

For an arbitrary SP spectrum the computation of the 
density of excited states of a fermionic system is a diffi- 
cult combinatorial problem for which no exact solution is 
known. We expect, however, that Eq. © gives correctly 
the contribution of the average part of the SP spectrum. 
What we are seeking here are the variations of the MB 
density due to fluctuations of the SP spectrum with re- 
spect to a perfectly ordered spectrum. A convenient way 
to express an approximate solution is by means of an 
inverse Laplace transform of Eq. J5J) . A saddle point ap- 
proximation of the resulting integrals yields 0,0 

p MB (A, E) = e s ^> T V 2^\V(p,T)\ , (4) 

valid in the degenerate gas approximation. The depen- 
dence on A and E in Eq. Q arises from the saddle point 
conditions that fix the value of the chemical potential p 
and temperature T of the gas 

Af(p,T)=A, £{p,T)=E, (5) 

The functions Af = —d£l/dp\ T and 

£ (p, T) = Q(p, T) + pM(p, T) + TS(p, T) (6) 

are the particle number and energy functions of the gas, 
respectively, with the entropy S(p,T) = —dVl/dT\^ and 
the grand potential Q = —T J dep(e) log [l + e^" e)/T ] , 
where pie) — Yli <K e — e ^ s the SP density of 
states. In terms of the previous functions, the 
determinant in Eq.QJ is defined as T)(p,T) — 
T 3 [dM/dp\ T d£/dT\„ - 8Af/dT\„ d£/dp\ T ]. All the 
necessary quantities involved in the computation of p MB 
are therefore defined from il(/i, T). 

The task consists in the evaluation of the entropy and 
the determinant at values of p and T that satisfy the 
conditions JSJl. As A and E are varied in Eqs.JSJ, in 
general p and T do not have a smooth and gentle be- 
havior, because the functions M and £ may have sudden 
changes due to the discrete nature of the SP spectrum. 
The difficulty lies in a proper treatment of the smooth 
part of the variations as well as the fluctuations with re- 
spect to it. To this purpose it is convenient to use the 
semiclassical approximation . The SP density of states 



2 

is decomposed into the sum of a smooth term (p, given 
by a Thomas Fermi approximation or Weyl series) plus 
oscillatory terms (p), p = ~p + p. Taking into account 
only the smooth part of the SP density of states, the 
calculation of p MB is relatively straightforward and leads 
to Bethe's formula. The oscillatory part depends on the 
primitive classical periodic orbits p (and their repetitions 
r), p — 2 J2 P Y^rLi Ap,r cos [rSp/H + v p>r ]. Each orbit is 
characterized by its action S p , stability amplitude A p ^ r , 
and Maslov index v p r . When inserted into the definition 
of the grand potential, and after integration with respect 
to the energy, the fluctuating part of O is given by |10| 

p r—l P 

Here t p is the period of the periodic orbit, and k[x) = 
x/sm\\{x) is a temperature factor that introduces the 
time scale t t = K/(ttT) conjugate to the temperature. 
This expression describes the departures of fl with re- 
spect to its mean behavior due to the fluctuations of the 
SP spectrum. By simple derivation of the grand potential 
the smooth and fluctuating part of all other thermody- 
namic functions can be obtained. Intensive functions like 
the chemical potential and the temperature should also 
be decomposed. We define their smooth parts by the 
implicit equations 

77(-p,T) = A, £(JI,T)=E, (8) 

whereas the fluctuating parts are p = p — ~p and T = 
T-T. 

Using Eqs.JSJ and the decomposition of the different 
functions, Eq.© may be rewritten as 

(T + f)S = E — H(p, T) - n(p, T)-(p + p)A . (9) 

If A is large, the smooth part of the grand potential can 
be expanded in power series of p and T. Keeping terms 
up to first order, using that in the degenerate gas approx- 
imation Q(p,T) ~ Q(/Z, 0) — (n 2 /6)p T , and neglecting 
variations of p, Eq.© becomes 

2 

(T + f)S = E-Tl{%Q) + ^-pT 2 +S{p,T) f- 

h(p,T)~pA. (10) 

This equation is valid for an arbitrary temperature. At 
zero temperature, the energy of the system is the ground 
state energy E , and Ea. fl7)| reduces to 

= E Q -Q(p, 0)-n( M o,0)- JIA . (II) 

~p remains constant if variations of p are neglected, but 
the fluctuating part still depends on temperature. As a 
consequence the chemical potentials are slightly differ- 
ent. The idea is to subtract Eq. Qllfl from I|1U|) . Before 
doing that we first need to properly analyze some of the 
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terms. The difference E — Eq = Q is the excitation en- 
ergy of the gas, as well as the term (7r 2 /6)p T = Q. 
The latter equation follows from the dependence of the 
energy on temperature in the degenerate gas approxima- 

2 

tion, £(/Z, T) ~ £(71, 0) + (tt 2 /6)p T , and the stationary 
phase condition (jSJ) for £ (and the corresponding one at 
zero temperature). It allows to express the temperature 
in terms of Q. Similarly, ~p is obtained by inversion of 
77 (ji, T) in (JSJ which, in the degenerate gas approxima- 
tion and neglecting variations of p, is independent of T . 

Returning to Eq.JTJJ, the term f(S - S(JI,T)) = TS 
is of second order, and is therefore neglected. After sub- 
traction of Ea. Hllf) and expressing T in terms of Q, the 
entropy may be expressed as 



S(A,Q) 



2 1 r 

r pQ+ f 



nfcO) -SlfcT) (12) 



-7T 2 p 



Q + n(-p,o) -n(-p,T) , (13) 



where, to lowest order, f2 has been evaluated at (ji, T) 
(a similar though improved accuracy is obtained using 
p, and T, cf. Ea.ljTo)! below). Equation (|T3l is valid if 

f2(/Z, 0) — f2(/Z, T) "C (J, and puts the result under the 
form of a backshifted Bethe formula. 

A more physical interpretation of Eas. (|12|) and (|13|l . 
that expresses the result in an entirely microcanonical 
language, is provided by the following connection. It has 
been shown previously [TTL Il2j | that, to leading order in 
an expansion in terms of ju, Q(/x,T) = £(-A, Q), where 
£(A,Q) are the (shell) fluctuations of the energy of the 
gas at a fixed number of particles and excitation energy. 
The shell corrections of Eq. (|13fl are therefore directly re- 
lated to the fluctuations of the energy of the system. 

A similar analysis can be performed for the determi- 
nant T>(fi,T) in Eq.Q. To leading order we find 



„ tt z 9 — 4 12 ^ 9 
V = —f T = —Q 2 



(14) 



recovering the same form obtained by Bethe (2iryV = 
\/48 Q, see Eq.(QJ). There exist oscillatory corrections to 
T>, but these are exponentially small compared to those 
associated with Ea. (|12f> . and we neglect them. 

Equation l|12|) is the central result of the paper. It ex- 
presses, together with Eq. l|14[) . the MB density of states 
in terms of the particle number A and excitation energy 
Q. It decomposes S into a smooth and a fluctuating 
part, S = S + S . The smooth part is given by the usual 
Bethe's result (corrections to it may be obtained by keep- 
ing higher order terms in p and T in the expansion). The 
novelty in Eq. (|12fl is the additional contribution of oscil- 
latory corrections, S = p(~p,Q) - ft(p,T)]/T. They de- 
scribe shell corrections to the MB density, ft depends on 
temperature only through the function n. The effect of 
this function is to exponentially suppress the contribution 
of periodic orbits whose period t p 3> r T |ToL Hl| . Since 



there is no suppression at T — (because r T — > 00), only 
orbits whose period r p > r T contribute to the difference 

Q(~p, 0) — T). At temperatures such that r r <C T m - m , 
where r min is the period of the shortest periodic orbit 
of the system, the term f2(p, T) becomes exponentially 
small, and only f2(ju, 0) remains. The shell correction 

S therefore decays as T . This decay contrasts with 
the more common exponential damping observed in other 
thermodynamic quantities |llj |. 

The behavior of fl (and therefore of S) strongly de- 
pends on whether /iorT are varied. As discussed above, 
a temperature variation modifies the prefactors of the 
summands in Q (through the function k), and therefore 
produces gentle variations of S. In contrast, A p . r , t p and 
S p depend on p. For large values of //, S p ^> h and 
the dominant variation with the particle number (or any 
other parameter that modifies the actions) comes from 
the argument of the cosine function in Vt. Rapid oscil- 
lations of S are therefore generically expected when the 
number of particles is varied. 

A clearer picture of how the fluctuations behave may 
be obtained through a statistical analysis. There are two 
relevant SP scales in the analysis, the mean spacing 5 
and the energy conjugate to the shortest periodic orbit, 
E c = h/r min . The ratio g = E c /S is typically much 
larger than 1 (for instance, g ~ A 2 ^ 3 in a three dimen- 
sional cavity). The most simple property of the fluctu- 
ations is (S) = 0, where the brackets denote an average 
over a suitable chemical potential window. This result 
is valid only to first order in the expansion; it can be 
shown that higher order terms contribute to a non-zero 
average. The next non-trivial statistical property is the 
variance (S 2 ), that may be computed using Eq.Q). The 
result is (S 2 ) = (1/2) J™ dx K (x , x K ) [I - k{x)] 2 / x i , 
where K(x, Xr) is the rescaled form factor of the SP 
spectrum (cf Eq.(36) in Ref.^l)- The latter function 
depends on the rescaled Heisenberg time xh = h~p/r T . It 
describes system-dependent features for x of the order of 
^min = T min/ T Ti while it is believed to be universal for 
x 3> x m - m . The universality class depends on the regular 
or chaotic nature of the dynamics, and on its symme- 
try properties. Taking into account the basic properties 
of K(x,xb), in chaotic systems we find three different 
regimes for {S 2 } as a function of T: 

(i) Low temperatures 2tt 2 T <C 6. In this regime (1S 2 ) = 
c a tt 2 T/5, where c 4 = 0.0609.... 

(ii) Intermediate temperatures 6 <C 27r 2 T <C E c . In this 
regime the size of the fluctuations saturate at a universal 
constant (S 2 ) = c 3 //3, where c 3 = 0.1023 ... and /3 = 1(2) 
for systems with (without) time-reversal invariance. 

(iii) High temperatures 27r 2 T ^> E c , The size of the 
fluctuations decreases with excitation energy, (S 2 ) = 
(n 2 (p,0))[l - 8 e- 27r2 T/ E -}/T 2 . After an exponential 
transient, a power-law decay (5 2 ) 1 / 2 oc T is obtained. 

The situation is different in integrable systems, where 



4 



25 



20 



IS 15 



A 



T=146.0T 8 = 1.45T C 
/\T=108.8T, = 1.08T e * ^ VA 



V 



T=48.7 T 5 = 0.48 T c 



1800 1850 1900 1950 2000 2050 2100 2150 2200 
A 

FIG. 1: The logarithm of the many-body density of states 
of A noninteracting fermions in a rectangular billiard of sides 

a — (I + VE)/2 and b = a -1 . Dots: numerical computation 
at three different excitation energies; solid curves: theoretical 
prediction 115H : dashed curves: smooth part Eq.©. 



only two regimes are found. At low temperatures the 
result is identical to that of chaotic systems. The differ- 
ence is that now the growth extends up to much higher 
temperatures, 2ir 2 T rs E c , without saturation. At that 
temperature the variance of the fluctuations is of order 
g. In integrable systems, the maximum amplitude of the 
fluctuations is therefore reached at 2ir 2 T ss E c , and its 
typical size is much larger than in chaotic systems. At 
high temperatures 2ir 2 T 3> E c the decay is almost identi- 
cal to that of chaotic systems (the coefficient 8 is replaced 
by a 12). 

As shown elsewhere T) is dominated, at any 

T, by the shortest classical periodic orbits. In contrast, 
the difference Q(7I, 0) — Q(~p, T) depends on orbits whose 



period t p ~t t . For temperatures 2ir 2 T <C E c the statis- 
tical properties of these orbits are universal, and corre- 
spondingly the probability distribution function of S is 
expected to be universal, in the sense that at a given tem- 
perature it should only depend on the nature of the un- 
derlying classical dynamics (regular or chaotic), and the 
symmetries of the system. This statement is supported 
by the fact that in the limit T^0,5« -80.(11, T) /dT. 
The probability distribution of the latter quantity was 
studied in Ref . JTT| i it was shown that it coincides at low 
temperatures with that obtained from a Poisson spec- 
trum for integrable systems and from a random matrix 
spectrum for chaotic ones. As the temperature is raised, 
the universality of the probability distribution of S will 
be lost for temperatures of the order or greater than E c , 
where system specific features are revealed. 

We have checked some of our predictions by a di- 
rect numerical counting of the MB density of states in 
a particular system. Figure 1 shows the results obtained 
for a gas of about 2000 fermions contained in a two- 
dimensional rectangular cavity, an integrable system. For 
each number of particles we compute the MB density of 
states at three different temperatures, measured in units 
of T s = 5/2n 2 and T c = E c /2tt 2 . The theoretical curve 
is computed according to the expression 

log(/Wp) = M/C/p) + [n(A*, 0) - n(fi,T)]/T , (15) 

where Q is given by the periodic orbits of the rectangle. 
Though the results are similar, u. instead of Jt is used 
to obtained a more accurate description of the numerical 
data. A systematic deviation is observed between theory 
and numerics if the corrections arising from the exact 
expression Eq.© are not included. Notice the extremely 
good accuracy of this equation, either for the average 
value of the density as well as for the fluctuations. 
The Laboratoire de Physique Theorique et Modeles 
Statistiques is an Unite de recherche de PUniversite Paris 
XI associee au CNRS. 
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